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Black Holes have always played a central role in investigations of quantum gravity. This in- 



cludes both conceptual issues such as the role of classical singularities and information loss, 



and technical ones to probe the consistency of candidate theories. Lacking a full theory of 
quantum gravity, such studies had long been restricted to black hole models which include 



some aspects of quantization. However, it is then not always clear whether the results 
$1, \ are consequences of quantum gravity per se or of the particular steps one had undertaken 
to bring the system into a treatable form. Over a little more than the last decade loop 
quantum gravity has emerged as a widely studied candidate for quantum gravity, where it 
^ ■ is now possible to introduce black hole models within a quantum theory of gravity. This 
makes it possible to use only quantum effects which are known to arise also in the full 
theory, but still work in a rather simple and physically interesting context of black holes. 
Recent developments have now led to the first physical results about non-rotating quantum 
black holes obtained in this way. Restricting to the interior inside the Schwarzschild hori- 
zon, the resulting quantum model is free of the classical singularity, which is a consequence 
of discrete quantum geometry taking over for the continuous classical space-time picture. 
This fact results in a change of paradigm concerning the information loss problem. The 
horizon itself can also be studied in the quantum theory by imposing horizon conditions 
at the level of states. Thereby one can illustrate the nature of horizon degrees of freedom 
and horizon fluctuations. All these developments allow us to study the quantum dynamics 
explicitly and in detail which provides a rich ground to test the consistency of the full 
theory. 
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1 Introduction 

Black holes in classical general relativity are, compared to other astrophysical objects, 
distinguished by the presence of singularities, where curvature and tidal forces diverge 
and where space-time stops, and horizons, which can separate off regions from causal 
contact from another region. Both properties have long been suspected to be changed in 
a quantum theory of gravity: Singularities denote points where the classical theory breaks 
down, and at least space-like ones which lie to the past or future of observers are supposed 
to be removed in a more complete quantum theory. Horizons, on the other hand, are still 
expected to play an important role also in quantum gravity. The horizon surface should 
at most be smeared out due to fluctuations in the causal structure on which the concept 
of horizons relies. For massive black holes (compared to the Planck mass) these horizon 
fluctuations should be negligible for most purposes such that the classical picture still 
applies. Instead of modifying the horizon on large scales, quantum gravity is expected to 
provide a microscopic picture which shows how to build a macroscopic horizon from Planck 
scale ingredients. If successful, this will then result in a statistical explanation of black 
hole entropy. 

In more detail, the main issues concerning black holes in quantum gravity are as follows: 

Singularities: Are they indeed removed and, if yes, what replaces them? There are 
arguments that not all singularities are equal, with space-like ones to be removed 
and time-like ones to persist in order to rule out unwanted (such as negative mass) 
solutions Also the issue of naked singularities and cosmic censorship arises in 
this context. 

Horizons: First of all, one has to see what an adequate definition of a horizon in quantum 
gravity could be. The original concept of the event horizon relies on the classical 
causal structure of all of space-time as well as the presence of singularities in the 
future. The quasi-local concept of isolated or dynamical horizons :2J uses much 
weaker assumptions about the structure of space-time such that it is better suited 
to a quantum treatment at least for large, semiclassical black holes which have only 
weak curvature at the horizon. For microscopic or primordial black holes, space-time 
even around the horizon cannot be treated as a smooth classical geometry with a 
classical causal structure. Here it is not clear if a quantum concept of horizon can 
even be applied. 

If there is an applicable notion of quantum horizon, the issue of black hole entropy 
can be analyzed. By identifying and counting quantum states uniquely characterizing 
a horizon of a given area one can compute black hole entropy and compare with the 
expected semiclassical Bekenstein-Hawking formula. Moreover, detailed pictures of 
the horizon structure and its fluctuations can be developed which shed more light on 
quantum gravity in general. If matter fields are present, the horizon should shrink 
from Hawking radiation which provides insights on how gravity interacts with matter 
at the quantum level. 
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1 INTRODUCTION 



Both: Systems such as black holes with singularities as well as horizons have led to much 
confusion in attempts to guess the outcome of quantum gravity from early glimpses 
obtained from mainly semiclassical considerations. This is most commonly expressed 
in the infamous information loss paradox according to which information falling into 
the singularity implies a non-unitary quantum evolution and thus presumably funda- 
mental limitations to knowledge [3]. These ideas obviously do not take into account 
what happens to singularities in quantum gravity and thus have to be revisited once 
a more complete treatment is available. 

All these issues probe different aspects of the full theory of quantum gravity and re- 
quire different techniques. A common feature, except for the entropy counting of isolated 
horizons, is that they are dynamical aspects such that the Hamiltonian constraint operator 
in a canonical quantization or an alternative evolution equation is essential. In particular, 
both black hole singularities as well as their horizons require inhomogeneous situations and 
an approximation by spatial homogeneity, which works well in cosmological cases, is not 
sufficient in general to grasp all the important physical aspects. This has the advantage of 
providing many non-trivial tests of quantum gravity which go beyond what is possible in 
homogeneous cosmological models. 

It certainly also implies that the treatment is more complicated, and indeed progress 
on the problems listed here has been mixed. The strongest results exist for the counting 
of black hole entropy of static or isolated horizons which has been derived in different 
approaches jU 13 E] • This has been possible since the isolation (or even extremality 
in allows one to ignore the complicated quantum dynamics and still compute the 
correct number of physical states. Moreover, only the horizon itself is important such 
that its inhomogeneous neighborhood does not have much influence. This changes if one 
also wants to study, e.g., horizon fluctuations since they are dynamical and require the 
neighborhood in which the horizon fluctuates. Thus, both the quantum dynamics and 
inhomogeneous configurations have to be handled, and there are not many results within 
a full candidate of quantum gravity so far. 

Similarly, the issue of singularities relies on dynamical aspects which for most of the time 
was too complicated to allow definitive conclusions as to whether or not singularities persist 
in quantum gravity. In the last few years, there has been progress on the homogeneous 
situation of cosmological singularities EH IH| which have been shown to be removed by 
quantum gravity • Analogous techniques are now also available for some inhomogeneous 
situations such as the spherically symmetric model f3J E] which is classically relevant for 
non-rotating black holes. This has led to an extension of the non-singularity statements 
from homogeneous models to the spherically symmetric one [T^]. Moreover, with new 
results about quantum horizons a consistent picture of quantum physics of black holes is 
emerging. 

This chapter is also intended as an introduction, by way of examples, to some of the 
techniques of quantum geometry with an emphasis on aspects which are typical for a loop 
quantization and essential for physical issues. The main theme will be the understanding 
of quantum dynamics in inhomogeneous situations and problems surrounding it. 
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2 Classical aspects of spherically symmetric systems 

A spherically symmetric metric is most easily written in polar coordinates (x, <p) and 
takes the form (with df2 2 = d$ 2 + sin 2 -My? 2 ) 

ds 2 = -N(x,t) 2 dt 2 + q xx (x,t)(dx + N x (x,t)dt) 2 + q w (x,t)dn 2 (1) 

where fields only depend on time t and the coordinate x of the 1-dimensional radial manifold 
B. This expression makes use of the lapse function N(x, t) and shift vector N x (x, t) which 
are prescribed by the slicing of space-time into spatial constant-t slices: coordinate time 
translations are generated by the vector field 

with the unit vector field n being normal to the slices. The spatial metric on those slices 
is then 

dg 2 = q xx (x, t)dx 2 + q w {x, t)dQ 2 (3) 

and extrinsic curvature 

K a b = \C n q ab , (4) 

which determines the conjugate ir ab = — ^y/det q(K ab — q ab K^) to the metric in a canonical 
formulation 16J, takes a similar form K = K xx (x,t)dx 2 + K w (x,t)dQ 2 . 

A well-known example is obtained by the spherically symmetric vacuum solution to 
Einstein's field equations, the Schwarzschild metric ^7] 

ds 2 = -(1 - 2M/x)dt 2 + - 1 dx 2 + x 2 dtt 2 (5) 

with the mass parameter M. It has the following properties: If we first restrict our attention 
to larger x > 2M, the metric is static since its coefficients do not depend on time and 
N x = 0. When x becomes large compared to the mass, i.e. if we approach the asymptotic 
regime far away from the black hole, the metric becomes asymptotically flat. The black hole 
region is characterized by the horizon which appears at x — 2M as a coordinate singularity 
in the Schwarzschild metric and can be defined in a coordinate independent manner as the 
outer boundary of a region where trapped surfaces, i.e. envelops of light rays which cannot 
expand outwards to infinity, occur. If we enter the black hole region through the horizon 
we notice that now t becomes a space-like coordinate since the tt component changes sign. 
The role of coordinate time is then played by x on which the metric coefficients depend. 
Thus, the interior is not static, but since the metric components now do not depend on 
the spatial coordinate t it is homogeneous (of Kantowski-Sachs form). 

2.1 Metric and triad 

The metric components q w = x 2 , q xx = (1 — 2M/x)~ 1 and ./V 2 = — g tt = 1 — 2M/x can be 
used to characterize the three different regimes of a massive black hole with mass M ^> 1: 



4 



2 SPHERICAL SYMMETRY 



At asymptotic infinity we have x ^> 2M ^> 1 and thus 

q w > i q X x ~ i • 

At the horizon we have x ~ 2M and 

q w > 1 q X x > 1 
while at the singularity we have ~ x <C 2M and 

|fe|<l N>1. 

In the latter relevant only if we approach the singularity on slices with t 

constant which are time-like inside the horizon. The lapse function, on the other hand, is 
the relevant metric component if we approach the singularity on slices which are space-like 
inside. 

These regimes of metric components can be used for a first glimpse on how a quantiza- 
tion may deal with the singularity or horizon. From cosmological models it is known that 
expressions for, e.g., curvature components can be modified when they become large, cut- 
ting off classical divergences (in isotropic cosmology they are all inverse powers of the scale 
factor [15] . or spin connection components in anisotropic models ^1]). Similarly here, some 
spin connection components contain information about intrinsic curvature. Their form can 
be obtained from the general expression (see, e.g., [12]) 
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where e\ are components of the co-triad (i.e. e l a e\ = q a b) and of its inverse. In spherical 
symmetry, co-triads take a special form just as the metric © does. Since it does not 
matter how a triad is rotated, it need not be exactly invariant under the rotation group 
acting on space, but it is enough for it to be invariant up to a gauge rotation. This is 
realized for co-triads of the form 



e l a Tidx a 



e x (x)r 3 dx + (ei(x)ri + e 2 (x)r 2 )d$ + (e 1 (x)r 2 - e 2 (x)r 1 ) smfidcp 
ex{x)r 3 dx + e (/3 (x)A(x)d^ + e ¥ ,(x)A(a;) sini?d<^ (7) 



where we use SU(2) generators Tj = — I Oj with Pauli matrices <jj, and A =: cos i]r 2 + sin t]T\ 
and A := exp(— fr 3 )A exp(^r 3 ) are defined to have unit norm in su(2), i.e. cosr/ = ei/e v 
and sin?] = —e-z/e^ with e^, = e\ + e\. Infinitesimal rotations of space now act by Lie 
derivatives on e with respect to superpositions of vector fields X = sin <pd# + cot cos (pd v , 
Y = —cos (fd# + cot $ sin (pd^, and Z = d v , while gauge rotations of the triad act by 
conjugation in su(2). We thus obtain explicitly 
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2.1 Metric and triad 



5 



showing that any rotation in space simply amounts to a gauge rotation of the triad. The 
corresponding metric is thus invariant under rotations, and indeed a co-triad ((7J) implies a 
metric of the form @ with 

_ 2 2,2 2 (o\ 

Qxx — £ x ' 1<P<P ~ e l ' e 2 ~ e ip ■ \°) 

A spherically symmetric spin connection takes the form 

T i r i dx a = T x r 3 dx + r v A r d$ + T [p A r sin Vdtp + r 3 cos tidip (9) 

where the last term must be added since a connection transforms differently from a co-triad 
under gauge transformations. Indeed, 

C x {r 3 cos "&d<p) = -r 3 
£y(t~3 cos$dy>) = — r 3 

and £z(t~3 cosi9d<£>) = such that we have the correct transformation of T with the same 
gauge rotation as above. 

The explicit formula (jBJ) applied to a spherically symmetric co-triad shows that 

T* = -rf , T v = -e' [p /e x , A r = A , (10) 

with A as defined for the co-triad, such that the (^-component is gauge invariant while T x 
is pure gauge. Modifications to classical behavior similar to those in cosmological models 
can now be expected, e.g., from the spin connection component = —^/q w '/ \fq^ when 
metric components become small. Classically, this expression diverges at small q xx , which 
can be changed in a quantum theory for the corresponding operator. Since, as we will see 
later, appears in the equations of motion, a modification here would change the behavior 
of solutions. Horizons of massive black holes would remain unmodified since there both 
metric components are large. The singularity, however, looks less clear: only for time-like 
slices does q xx become small, indicating a modification and the possibility of removal of 
the singularity. But if we approach the singularity on space-like slices with x constant, 
in the interior A^ 2 (playing then the role of q xx ) remains large which does not suggest 
modifications. Indeed, the slices then are homogeneous and vanishes identically which 
means that we will need another measure for the removal of singularities in this case. 

At asymptotic infinity, however, we would encounter severe problems since q xx is close 
to one at which point modifications can already be noticeable, spoiling the classical limit 
of the theory. This is a sign of using the wrong variables since the modification is a 
consequence of quantum effects, and the success of a quantization can depend significantly 
on the choice of fundamental variables. Indeed, there are variables better suited to a 
demarkation of the different regimes than the metric. This is in particular the case for 
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2 SPHERICAL SYMMETRY 



the densitized triad defined by Ef = ef | det(e^)| where ef is the inverse of the co-triad e 
compatible with the metric. A spherically symmetric densitized triad is of the general form 

E = E x (x)r 3 sin#— + {E\x)t 1 + E 2 {x)r 2 ) sintf- + {E\x)t 2 - E 2 (x)r 1 ) — (11) 

written down as an su(2) valued densitized vector field. The gauge invariant components 
are E x and (E^) 2 = (E 1 ) 2 + (E 2 ) 2 whose relation with the metric components is 

\E X \ = q w E* = Jq~q^ (12) 

(note that E x can be positive or negative depending on the orientation sgn det E = 
sgn E x (E lfi ) 2 of the triad). The angular components have the same internal directions 
A and A as the co-triad. 



For the Schwarzschild solution with \E X \ = x 2 and E^ = xj a/1 — 2M/x we now have 
the following behavior: At asymptotic infinity 

|£ x |»l 

at the horizon 

\E X \ > 1 E^ > 1 

and at the singularity 

|-E X |<1 £^<1. 

Thus, irrespective of the approach to the singularity, the behavior is just as needed for 
unmodified classical behavior far away from the black hole all the way up to the horizon, 
while inverse triad components, such as the spin connection component 

= -(E X )'/2E V (13) 

will be modified at the singularity with small E v . 



2.2 Basic variables 

For detecting the classical singularity it seems much more reliable to use the densitized 
triad rather than the metric, which is also the case in homogeneous models with an explicit 
impact on the removal of singularities ^U]- Indeed, the densitized triad as a basic variable 
is important in other ways, too: it arises naturally when one attempts to quantize gravity 
in a background independent manner. These two issues, the fate of classical singularities 
and background independence, are superficially quite different but turn out to be deeply 
related. 

Most recent progress in a background independent quantization of general relativity 
has come after a reformulation in terms of Ashtekar variables |2U| l2Tj where the densitized 
triad Ef plays the role of a momentum canonically conjugate to the Ashtekar connection 
A l a = F l a — r yK l a with the spin connection T l a as a function of Ef via (JBJ), extrinsic curvature 



2.2 Basic variables 
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K l a = e\Kab and the Barbero-Immirzi parameter 7 > [221- The extrinsic curvature 
components here make A l a canonically conjugate to Ef, while the spin connection provides 
A l a with the transformation properties of a connection. This reformulation thus casts 
general relativity as a gauge theory and does not only bring it formally closer to other 
interactions but also leads to a direct way for a background independent quantization. 

Usually, a field theory would be quantized by smearing the fields with test functions 
over 3-dimensional regions so as to make their classical Poisson *-algebra well defined. 
For instance, a scalar <fi with Lagrangian y/det q{\4> 2 + \q ab d a (j)db(t> + V{4>)) on a back- 
ground metric q a b (assuming lapse function N = 1 and shift vector N a = 0) has mo- 
mentum = \J det q<p which transforms as a density (which is often ignored when the 
background metric is fixed as, e.g., Minkowski space). This has the singular Poisson re- 
lations {4>(x) , Ptj>(y)} = S(x,y). However, if we smear the fields with test functions / and 
g on space to obtain <p[f] := J y/det qf(x)(p(x)d 3 x and p<j,[g\ '■— f g(x)p ( f,(x)d 3 x we ob- 
tain the well-defined Poisson algebra {</>[/],£></,[#]} = / y/detqf(x)g(x)d 3 x. This does not 
contain 5-functions, but does depend on the background metric q which is not available 
for a background independent formulation of gravity. The very first step of a background 
independent quantization of general relativity, therefore, has to face the problem that the 
physical fields, with the metric or densitized triad among them, need to be smeared for a 
well-defined algebra to be represented on a Hilbert space, but that a background metric 
must not be introduced. 

For a scalar, there is a simple way out: as is easily verified, we still obtain a well- 
defined algebra if we only smear p^ for which we do not need a background metric since it 
is already a density. Similarly, in the case of gravity we can evade the problem in Ashtekar's 
formulation since with connections and densitized vector fields as basic variables there is a 
natural, background independent smearing leading to a well-defined algebra: Instead of 3- 
dimensional smearings for all basic fields we use a 1-dimensional smearing of the connection 
and a 2-dimensional one for the densitized triad, giving rise to holonomies 



through surfaces S. (We use the tangent vector e a to the curve e and the co-normal n a to 
the surface S, both of which are defined without reference to a metric.) 

It turns out that this smearing is sufficient for a well-defined classical Poisson algebra 
which even has a unique diffeomorphism invariant representation [2B1 I2H 123 1213 12Z| • This 
representation defines the basic framework of loop quantum gravity [2E1 12H1 CUB EI] ■ States 
are represented usually in the connection representation ip[A] on which holonomies act as 
multiplication operators and fluxes as derivative operators. This can all be done rigorously 
thanks to a rich structure on the infinite dimensional space of connections which is under 




(14) 



along edges e in space, and fluxes 




(15) 
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2 SPHERICAL SYMMETRY 



much better control than the space of metrics. As a consequence, flux operators have 
discrete spectra implying a discrete structure of spatial geometry [321 HUH Ell which is also 
realized in symmetric models |85}l36j. Moreover, since flux spectra are discrete and contain 
zero, there are no densely defined inverse operators. Instead there are techniques |HZ| which 
allow one to quantize co-triad or other inverse components of the basic E? by operators 
which reduce to the inverse in a classical regime but modify the classical divergence at small 
values. This has already been described and used above for the spherically symmetric 
spin connection component. Here, it is important that those expressions are taken as 
functions of the densitized triad components and not metric components. These effects 
come from properties of flux operators as basic operators in a background independent 
formulation which relies on the densitized triad as basic variable and so far is not known 
in a metric formulation. Indeed, as observed before, the densitized triad is much better 
suited to separate the classical singularity from other regimes such that modifications are 
only expected there. 

2.3 Dynamics 

Up until now we have discussed kinematical properties of the spherically symmetric system. 
The dynamical behavior of triad and connection (or extrinsic curvature) components is 
dictated by the Hamiltonian constraint 

H[N] = (2G)- 1 [ dxN(x)\E x \-^ 2 ((K 2 E* + 2K V K X E X ) + (1 - T 2 )E^ + 2V ' E x ) (16) 
Jb 

in terms of the spin connection component as before and the extrinsic curvature com- 
ponents in 

K = K x (x)r 3 dx + {K x (x)t x + K 2 (x)r 2 )d$ + (K^xfo - K 2 (x)n) sin&dip (17) 

where again only K x and K 2 = K\ + K\ are gauge invariant. In addition, there is the 
diffeomorphism constraint 

D[N X ] = (2G)- 1 [ iV x (x)(-2£^ + K X E X ') . (18) 
Jb 

Physical fields (K x , E x ; K^, E v ) have to solve the constraint equations H[N] = = D[N X ] 
for all functions N and N x on B (except for possible boundary conditions which we ignore 
here) and evolve in coordinate time according to Hamiltonian equations of motion E x = 
{E x , H[N] + -DfiV^']}, etc. to be computed with the Poisson relations {K x (xi), E x {x2)} = 
—2GS(xi,x 2 ), {K v {x\)i E lf (x2) = —Gd(xi,X2)- For the triad components this gives 

E x = 2NK v y/\EF\+ N X E XI (19) 
E v = N(K V E V + K X E X )\E X \- 1/2 + (N x E ip y (20) 

which, when solved for the extrinsic curvature components, agrees with their geometrical 
definition via 

K = e\K ah = (2Ny l e\C dt „ N * dx ei a e{ (21) 



2.3 Dynamics 
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from (J3J) and (0). Evaluating this for a spherically symmetric co-triad (J7|) or densitized 
triad (fTTj) indeed gives spherically symmetric components 

K x = N-\e x -(N x e x )') , = N~\e v - N x e' v ) 

and the same internal directions A K = A, A K = A as those of the triad. The extrinsic 
curvature components then have Hamiltonian equations of motion 

K x = -NK V K X \E X \~ 1/2 + \NKlE v \E x \-* /2 + {N X K X )' 

+\N\E X \- 1/2 [E^E*^ 1 - i(E x ') 2 (\E x \E' p )- 1 - E x ' E*\E^y 2 + E^'iE^y 1 ) 
+ ±N' (V'^r^T 172 - 2^/\E^\E v '{E lp y 2 ^ + N" ^/IeFKE^)- 1 (22) 

K v = -INK^EX 1 ' 2 + N x K'^ 

+ \N\E x \- l ' 2 (\(E x ') 2 (E^)- 2 - 1) + \N' ^\eF\E x \E*)- 2 . (23) 

These coupled non-linear equations are difficult to solve in general, but the Schwarzschild 
solution can easily be reproduced by assuming staticity: K x = = N x = which al- 
ready implies that the diffeomorphism constraint is satisfied. Equations (}2*2*|) and (}2*3*j) then 
assume the form of consistency conditions for the lapse function in order to ensure the 
existence of a static slicing. Both conditions are identically satisfied for a lapse function 

N oc E x '/E v (24) 

using that E v and E x are subject to the constraint equation 

(rj - l)E v - 2T' lp E x = (\(E X ') 2 / (E*) 2 - 1)^ + [E x 'jE^)'E x = 

following from (fTBj) with K x = = K v and from f)13|) . 

It remains to solve this constraint for E x and . If we choose our radial coordinate 
such that \E X \ = x 2 , this simplifies to a differential equation 

-2x 3 E lp ' + 3x 2 E v - (E^f = 

whose solution E v {x) = x(l + c/x)~ l l 2 is the Schwarzschild component for E v with c = 
— 2M, which then also reproduces the correct lapse function from (|2~Uh 

This shows how the dynamical equations appear in a canonical formalism, and also how 
special the simplicity of the static Schwarzschild solution is. With slight modifications to 
the equations, e.g. coming from quantum modifications, the assumption of staticity will no 
longer be consistent since two conditions K x = = K v have to be satisfied by only one 
function N. Thus, quantum corrections are expected to change the static behavior of the 
classical solution, even though it would come from only small changes outside the horizons 
of massive black holes. What this means for the inside where quantum effects dominate 
around the singularity has to be analyzed by direct methods from quantum gravity. 
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3 QUANTIZATION: OVERVIEW 



3 Quantization: Overview 

Even though the vacuum spherically symmetric system has only a finite number of physical 
degrees of freedom given by the black hole ADM mass and its conjugate momentum jHHl 
EHIED], a Dirac quantization requires field theory techniques in order to deal with infinitely 
many kinematical degrees of freedom. Almost all of these degrees of freedom will then be 
removed by the Hamiltonian constraint which acts as a functional differential or difference 
operator. Thus, many of the field theoretic aspects of the full theory can be probed here 
which also implies a corresponding level of complication. So far, the system is not fully 
understood in a loop quantization even in the vacuum case, and other techniques which 
can be applied more easily to this system do not allow definitive conclusions about the 
singularity. It is therefore necessary at this stage to refer to approximation methods. 
These methods allow different glimpses which one can then try to bring together for a 
consistent picture. Here, we briefly collect different classes of approximations, which will 
be described in more detail in the following sections. 

3.1 Homogeneous techniques 

Currently, loop techniques for homogeneous geometries, following techniques introduced in 
[JT], are fully developed to a degree that one can analyze properties of physical solutions. 
(The main open issue is the physical inner product, about which not much is known even in 
the simplest cases jJ2JH31-) There are explicit expressions for the most important operators 
such as the volume operator [33], matter Hamiltonians or the Hamiltonian constraint [3 
ITU] ITT] which is a big advantage compared to the full theory where the ubiquitous volume 
operator cannot be diagonalized even in principle. The constraint equation takes the form 
of a difference equation for the wave function in the triad representation which explicitly 
shows how one can evolve through the classical singularity. Moreover, one can define 
effective classical equations with diverse correction terms j33J H5] Ell E3 EHl HH ED] • They 
capture the main quantum effects jUJ EJJ E2 EI3 EH EI3 EHl EH EH] and can be analyzed 
more easily than the quantum difference equation directly (see, e.g., |qTH IfiU] IdT] lo^] Id^ IM] 

ESI)- 

In some cases these effective classical equations provide an intuitive explanation for the 
removal of singularities since they display bouncing behavior of a cosmological solution. 
This can also be used to model the case of matter collapsing into a black hole. As a model, 
the ball of matter can be assumed to be homogeneous such that the collapse of the outer 
shell radius is described by effective equations for an isotropic system. These equations 
are modified at small scales, i.e. when the ball collapses to a certain size. In the modified 
regime there are matter systems which show a bounce, which now can be interpreted as the 
collapsing matter parts repelling each other and bouncing back after maximal contraction. 
So far, this is not much different from a bouncing universe and indeed described by the 
same equations. The difference is that the matter ball does not present the full system, but 
that there is also the outside. Without specifying the matter content there, one can try 
to match the interior to a generalized Vaidya metric outside allowing for matter radiated 
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away. This allows to study the formation or disappearance of horizons which may or may 
not shield the bounce replacing the classical singularity [66 . 

Limitations of these techniques are that only the interior carries quantum effects, while 
the outside is described by a generalized Vaidya metric of general relativity. Some quantum 
effects are transported to the outside by matching to the effective interior, which then enter 
the Vaidya solution effectively through a non-standard energy momentum tensor. This still 
shows possible changes in the behavior of horizons, but is of course more indirect than a 
complete inhomogeneous analysis. 

A different approach using homogeneous techniques only applies to the Schwarzschild 
solution which is homogeneous inside the horizon. One can then describe the interior by 
a quantum equation which as in cosmological CcLSGS, IS cl difference equation not breaking 
down at the classical singularity. Also here we thus obtain a mechanism to evolve through 
a classical singularity, and there are many more non-trivial aspects which only arise in a 
loop quantization and show its consistency [HZ]. In particular, the singularity is removed, 
but the horizon which presents another boundary to the classical interior remains. 

3.2 Extrapolation 

The previous analysis provides a picture of a non-zero Schwarzschild black hole interior 
which one can now extrapolate in two ways: The non-singular interior first has to be 
embedded in a full space-time which can happen in several different ways. Moreover, for 
a realistic black hole this must be generalized to the presence of matter. While there are 
many gaps to be filled in by detailed constructions and calculations, one can already see 
different implications for the issue of information loss [68 . 




3.3 Inhomogeneous techniques 

Operators for the spherically symmetric system (with or without matter) are now available 
explicitly at a level similar to that in homogeneous models [T3| I14j. In particular, there is 
a similar simplification in the volume operator which translates to matrix elements of the 
Hamiltonian constraint also being known explicitly. However, the constraint is much more 
difficult to analyze since it now presents a functional difference equation in infinitely many 
kinematical variables. The construction and regularization of the constraint is more subtle 
compared to homogeneous cases, but similar to the full theory where there are different 
versions. These can then be studied explicitly and their physical implications analyzed, 
leading possibly to conclusions as to which operator is most suited for the full theory. 

Even though the singularity issue is not yet solved in generality, there are indications 
that a mechanism similar to that in homogeneous models is at work. This will then 
provide a large class of systems where one and the same mechanism, derived from basic 
loop properties, provides a removal of singularities in non-trivial ways. 

There are regimes where the constraint operator can be approximated by a simpler 
expression. Interestingly, this is true in particular in the neighborhood of isolated or slowly 
evolving horizons such that horizon properties such as fluctuations and its growth 




12 



4 HOMOGENEOUS TECHNIQUES 



from infalling matter or shrinking from Hawking radiation can be analyzed. Moreover, the 
regime provides perturbation techniques which allow us to study general properties of the 
constraint operator and matter Hamiltonians. 

3.4 Full theory 

The full theory has a rigorous quantum representation j7_Q] and well-defined candidates for 
the Hamiltonian constraint [7j]. Understanding the dynamics in general is certainly very 
complicated, and even computing matrix elements of the constraint is involved. Most full 
results which contribute to the physical picture are thus non-dynamical: Spatial geometry 
is discrete [321 E31 EH] as a characteristic of the full quantum representation, and there are 
well-defined quantum matter Hamiltonians [37]. Black hole (and cosmological) horizons 
can be introduced as a boundary provided they are isolated [2j. This condition ensures 
that the dynamics at the boundary is not essential and allows the correct counting of black 
hole entropy [7J [H] • 

4 Homogeneous techniques 

In the Schwarzschild interior r < 2M one can choose a homogeneous slicing such that the 
metric is of the Kantowski-Sachs form 



with T = r, R = t and a lapse function N(T) 2 = T/(2M — T). The spatial metric is 
then related to a homogeneous triad of the form (jll|) where E x and E v are constants 
on spatial slices. Their conjugates are given by Ashtekar connection components of the 
general spherically symmetric and homogeneous form 



which in this case are simply proportional to extrinsic curvature components K x = —A x /^f 
and K v = —A ip / r y since V = T 3 cos"&dif from Q and (jl(J|) with homogeneity. Moreover, 
A A = A (as defined for the triad) follows from the Gauss constraint. Since A is constant 
in a homogeneous model and subject to gauge rotations, we will fix it to A = t-i in this 
section, such that A = t\. The symplectic structure for the 4-dimensional phase space is 
determined by {K x , E x } = -2G, {K V ,E^} = -G. 

4.1 Quantum representation 

Loop quantum gravity is based on spin network states which are generated by holonomies 
as multiplication operators. Similarly, homogeneous models in loop quantum cosmology are 
based on a representation [72] which emerges from holonomies of homogeneous connections 
and which turns out to be inequivalent to the usual Schrodinger representation used in a 



ds 2 = -N(T) 2 dT 2 + (2M/T - l)dR 2 + T 2 dtt 2 



(25) 



A i a r i dx a = A x r 3 dx + A lf A A d$ + A^A A sin^dy? + r 3 cos$d<^ 



(26) 
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Wheeler-DeWitt like quantization. For the Kantowski-Sachs model an orthonormal basis 
of states is given by the family 

(K v ,K x \fx,u) =exrj(-^(fxK v + uK x )) /i,z/eM,/i>0 (27) 

such that the kinematical Hilbert space is non-separable. (There are arguments to reduce 
this to a separable Hilbert space as in [73] using properties of observables [H].) One can see 
one of the basic loop properties that only exponentials of connection or extrinsic curvature 
components are represented directly, but not the components themselves: It is clear that, 
e.g., exp(— i r yK x /2) acts directly shift operator 

exp(— i-fKK x )\fi, v) = v + 2k) (28) 

but since this operator family is not represented continuously, this does not allow us to 
obtain an operator for K x by differentiation. Indeed, 

(ji, u\exp(-ijKK x )\fi, u) = (n, u\fi, v + 2k) = 5 , K 

is not continuous at k = 0. This is different from a Wheeler-DeWitt quantization where 
extrinsic curvature components would be basic operators represented directly. Instead, 
here we have to express those components through holonomies such as exp^R^) = 
cos{\^K x ) + 2t3 sia(^jK x ) and use the action 

cos(i T iQ|/^) = i(|/i,z/ + l) + |/i,z/-l» (29) 
sin(§ 7 iQ|/^) = ^ - 1> - 1^, zy + 1>) . (30) 

Another difference to the Wheeler-DeWitt representation arises for triad operators 
which in the Wheeler-DeWitt case would be simply multiplication operators on a wave 
function in the metric representation and thus have continuous spectra. In the loop case, 
however, the triad operators 

P 2 r) * P 2 ri 

EX = i ^pd fr, = i h0 (31) 



with the Planck length £p = \Z87rGh have the basis states (j27J as eigenstates 

E x \u., v) = ±ry%v\», v) E*>\n, v) = ^j£ 2 Mfi, v) (32) 

and thus discrete spectra (i.e., normalizable eigenstates). Again, this is different from 
the Wheeler-DeWitt quantization but directly analogous to full loop quantum gravity. In 
particular the volume V = AnE^ \J\E X \ has a quantization with discrete spectrum with 
eigenvalues 

= 27r( 7 4/87r) 3 /V VW\ ■ (33) 
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4.2 Inverse triad components 

It is often necessary to quantize inverse powers of densitized triad components, for instance 
for matter Hamiltonians or curvature components. Since the basic triad operators have 
discrete spectra containing zero, they do not have densely defined inverses which could 
otherwise be used for this purpose. Nevertheless one can proceed, and in the end have 
regular properties, by rewriting the classical inverse in an equivalent way and quantizing the 
new expression |37|. We demonstrate this for the spatial curvature given by 3 R = 2/\E x \ 
for which we need an inverse of E x . This can be taken as a measure for the classical 
singularity where it diverges. Since there is no direct way of quantizing this expression via 
an inverse of E x we first write 

where the first step replaces the inverse power of E x by only positive powers occurring 
in V at the expense of introducing K x for which we do not have a direct quantization. 
Nevertheless, in the next step we obtain an equivalent expression which only contains 
exponentials of K x which we can quantize directly. Using the volume operator and turning 
the Poisson bracket into a commutator then yields a densely defined operator 

El^£l = _ZL^trr 3 e-^[e^,t>] (34) 
Ai 



i v sin(±7fQV cos(\^K x ) - cos (±7^)^ sin (±7^)) 



with eigenvalues 



(V,, v+ i - V^-i) = ky^-Ky/W + M - VW^M) • (35) 



Since E^ has eigenvalues 7^/1/ 167r, we can write 
sgn E' 



fi, v) = ( 7 4/8vr)- 1/2 ( VW+T\ - VW^M) (36) 



which has the expected behavior for \u\ ^> 1 but behaves very differently from the classical 
expectation for small \u\. 

Taking this as a measure for the singularity indicates that it is removed in quantum 
gravity since the eigenvalues remain finite even when v — at the classical singularity. 
Nevertheless, a final confirmation of an absence of singularities can only come from consid- 
erations of the dynamics which must allow us to evolve further even when we reach a point 
corresponding to the classical singularity. Only then can we conclude that the singularity 
boundary of space-time has been removed. 
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4.3 Dynamics 

The spherically symmetric Hamiltonian constraint (J 16)) can be used to find the expression 
for the homogeneous Kantowski-Sachs interior 

H[N] = (2G)- 1 N\E X \- 1/2 {{Kl + 1)E^ + 2K V K X E X )) (37) 

where we used = with homogeneity in ()13j) . There are different terms in this ex- 
pression, those quadratic in K and the i^-independent one which comes from the spin 
connection in the curvature of the Ashtekar connection. In the full theory there would 
only be curvature components of A\ in the Euclidean part e^F^E^E^ of the constraint, 
which can be represented via holonomies by using 

^F4(x)r l - = (/i a -l)/A + 0(A) 

where a is a small loop of coordinate area A and with tangent vectors si and S2- For small 
A in a limit removing a regulator one can use h a as an excellent approximation for the 
curvature components, and stick this together with quantizations of the triad components 
to obtain a quantization of the constraint [TTj. This is different in a homogeneous context 
(or in any symmetric model where some directions are homogeneous) because we have 
only exponentials of connection components, but not holonomies with an adjustable edge 
length that shrinks in a continuum limit. Nevertheless, since the constraint operator in 
the full theory is based on holonomies quantizing the F-components, this has to be the 
case also for symmetric models related to the full theory. The only possibility to use h a as 
a good approximation is then given when the arguments of the exponentials are small in 
semiclassical regimes where the classical constraint is to be reproduced. In other regimes, 
one does not expect the classical constraint to be of any value for guidance and in fact 
usually obtains strong quantum corrections. 

In a semiclassical regime one has small curvature such that the extrinsic curvature 
components can be assumed to be small when checking the classical limit of the constraint. 
However, Ashtekar connection components are not necessarily small since for them also the 
spin connection plays a role. Here, another difference to the full theory arises: while in 
general spin connection components do not have coordinate independent meaning and in 
fact can be made arbitrarily small in any neighborhood, some of the components (such 
as (|13|) in the spherically symmetric model) obtain invariant meaning in a symmetric 
context where only transformations respecting the symmetry are allowed. Usually, unless 
the model has flat symmetry orbits such that the spin connection vanishes, one cannot 
expect the components to be small even in semiclassical regimes. This requires a special 
treatment of the spin connection in symmetric models, which is possible in a general manner 
P"H IT%j . For this reason we have started the quantization in this model with extrinsic 
curvature components and will also use them instead of Ashtekar connection components 
in holonomies when constructing the Hamiltonian constraint. 

One may ask what the relation to the full theory then is where holonomies of the 
Ashtekar connection are basic, while holonomies of a tensor such as extrinsic curvature 
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cannot even be denned. The arguments presented before explain why extrinsic curvature 
is important to analyze the classical limit, but this does not show the contact to the 
full theory. This will be much clearer in inhomogeneous models which are in between 
homogeneous ones and the full theory. Here, we will have directions along symmetry 
orbits, for which the techniques just described will apply, and inhomogeneous directions 
for which we will use holonomies of the Ashtekar connection as in the full theory. As we will 
discuss in more detail in the quantization of the spherically symmetric model, all this fits 
into a general scheme which allows to derive expressions in all different classes of models. 
We can now express the terms quadratic in curvature components via holonomies, such 

as 

2 

and 



Kl + 1 = — tr r 3 (e- i ^ 71 e-^ 75 e^ 71 e^ 75 + 7 2 5 2 r 3 ) + 0(5) 



K X K V = J^t ITl (e- s ~i KxT2 e- & ~i K ^e &lK * T2 e &lK « TZ ) + 0(5) . 
r y 2 5 I 

Triad components, together with Pauli matrices in the traces, can be obtained in the right 
combinations from the Poisson brackets 

p"P 1 
■7- — r -6~/K x T 3 f r SjK x r 3 vi 



as already used for (J3JJ), and 



47170 G 



In all expressions, besides the volume V = An J dx^\E x \E' p only holonomies h& '■ = 
e -7&K x r^ h (S) ._ and h (S) ._ e -i8K v n of the S y mme t r ic Ashtekar connection (jUJ), 

expressed through extrinsic curvature components, occur which can be quantized directly. 
In a more symmetric form, which as we will see later also applies in general, we write 

(K 2 n + l)E* + 2K X K W E X 1 ,,00 

' ' .tr((/*VSV +l 2 5 2 r 3 )h x {h-\V} 



+h x h#h~ 1 h# 1 h v {h~ 1 , V} + h^hji^h^h^h^ 1 , V}) 
£ e IJK tr ((hjhjhfhj 1 - 7 2 5 2 F(r) /J )/ ix {/ i - 1 , V}) 



' UK 

with the curvature components F(T)u of the spin connection, i.e. here 

F(T) = dr = -r 3 sin tfdtf A dip 

such that only F(T) &!fi := i Xlfi ix^F{T) = — t 3 appears, with the symmetry generators 
X$ = d$ and X v = (sin??)" 1 ^. 
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Quantizing and evaluating the action explicitly through the action of basic operators 
leads to a constraint operator of the form 

= -iW( 7 3 * 3 G4)- 1 ^e ZJJr tr [ihfhf hf^hf' 1 - y^F^jj^ [h^~\ V] 

UK 

3*W2 w f Q . 5 1 k <p 8<yK v 5yK x ____ 5yK x 



-2iN('f&Gli>) sin cos sin cos 

. <h^ T > fry-K*. hK v y . SyK v 
x sm 21 V cos 21 — cos V sm - 



2 (57^ 2 (57^ 2 r2 

+ ( 4 sm — — — cos — — - + 7 d 

, . 5jK x y 5yK x 5jK x y . 5yK x . . 
x I sm — - — 1/ cos — cos — - — V sm — - — (St 



where 5 > is regarded as a parameter analogous to the edge length in the full theory. 
From the holonomy operators one obtains shifts in the labels when acting on a state 
\fi, v) which in the triad representation given by the coefficients t/> M) „ in a decomposition 
1^) = v ^n,v\^i v) leads to the difference equation 

= (H (5) ij)^ = 25y/\v + 25\(lP tl +28,v+2S - ^-2S,U+2S) (39) 

mVW + &\ - y/\v-6\) ((/i + 45)^+45,, - 2(1 + y 2 5 2 )fi^^ u + 45)^_ 45 ,,) 

-25yJ\v - 251(^+25,^-25 - VV-25,z;-25) • 

We are now in a position to analyze whether or not there is a singularity in the quantum 
theory. There are a few key differences to the usual classical formulation, the first one 
coming from the fact that we are using triad variables. Compared to a metric formulation, 
this provides us with an additional sign factor sgn E x determining the orientation of space. 
Accordingly, there are two regions of minisuperspace separated by the line E x = where the 
classical singularity would be. We have already seen that the classical divergence of inverse 
powers of E x does not occur in a loop quantization, but the real test of a singularity can 
only come from the dynamics: starting with initial values in one region of minisuperspace 
we need to find out whether we can uniquely evolve to the other side, through the classical 
singularity. In the quantum theory this is done for the wave function which we can prescribe 
for sufficiently many initial values at some v > and additional boundary values at \x = 
so as to provide a good initial value formulation for the difference equation as described 
in detail in ^21- One can then see by direct inspection that indeed this will uniquely fix 
the wave function not just at positive v where we started, but also at negative v, at the 
other side of the classical singularity. Thus, quantum geometry automatically allows us to 
evolve through the classical singularity which therefore is removed from quantum gravity. 
Intuitively, the region of negative v corresponds to a region of a space-time diagram at the 
other side of the singularity, as sketched in Fig. ^ which therefore is no longer a boundary 
but a region of high curvature where the classical theory and its smooth space-time picture 
break down [UTj . 
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Figure 1: Interior of a Schwarzschild black 
hole with the quantum region (hatched) re- 
placing the classical singularity. This allows 
to extend space-time to the new upper re- 
gion. How these regions are embedded in a 
full space-time is left unspecified here. 



There are many basic aspects which are playing together in just the right way for this 
result to hold true. They all come directly from the loop quantization and are not put 
in by hand; in fact, they had been recognized as essential for a background independent 
quantization a few years before their role in removing classical singularities emerged. The 
loop representation is important in two ways since via discrete triad spectra it leads to 
the kinematical results of non-diverging inverse powers of densitized triad components, 
and through the representation of holonomies to the dynamical constraint as a difference 
operator. Moreover, the theory is based on densitized triads which, as discussed before, has 
consequences for the position of classical singularities in minisuperspace important for how 
one can evolve through them. This automatically provides us with the sign factor from 
orientation and thus a region beyond the classical singularity. Still, also the dynamical 
law has to be of the right form for an evolution to this other side of minisuperspace to be 
possible. 

Thus, we have a few essential effects which automatically come from a loop quantization. 
Once recognized and identified, they can easily be copied in other quantization schemes 
inspired by loop quantum gravity and cosmology. However, in such a case one has to guess 
anew in each model what the relevant basic properties would be since there is no underlying 
scheme for guidance. With the loop quantization we have such a general scheme which 
just needs to be evaluated in different models. Only then can the results be regarded as 
reliable expectations for quantum gravity, rather than possibly artificial consequences of 
choices made. The sign of triad components, for instance, appears automatically and then 
gives rise to the additional side of the classical singularity to which we can evolve. In loop 
inspired approaches without a link to the full theory, however, the sign is introduced by 
hand by extending the range of metric variables to negative values. While this leads to 
similar results in isotropic models [TSj, except that the geometrical meaning of the sign 
remains unclear, there are differences in the black hole interior [TH]. In particular, this 
approach would suggest that even the horizon can be penetrated by the homogeneous 
quantum evolution despite the fact that space-time becomes inhomogeneous outside. This 
problem does not occur in the quantization described here since the horizon remains a 
boundary (corresponding to \i = 0). 
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4.4 Effective dynamics 

The non-singular quantum dynamics is obviously very different from the classical one even 
though they can be shown to agree in classical regions at large densitized triad components 
and small curvature |77j. In between, there is a regime where equations of motion of the 
classical type, i.e. ordinary differential equations in coordinate time, should be able to 
describe the system even though quantum effects are already at work. One can think of 
these equations as describing the position of wave packets which spread only slightly in 
semiclassical regimes jlHl EEJ HZ] • Quantum effects will then provide modifications, e.g. 
where inverse powers of densitized triad components occur in a matter Hamiltonian which 
are replaced by regular expressions in quantum geometry. This provides different means 
to calculate implications of quantum effects which can so far be done in homogeneous 
situations. 

For instance if we assume the distribution of a matter system collapsing into a black 
hole to be isotropic, its outer radius Ra(t) is described by a solution a(t) to the Friedmann 
equation, with R being the coordinate radius where we cut off spatial slices from a closed 
FRW model. If we choose a scalar </> with potential V(4>) and momentum p^, we have the 
Friedmann equation 

a(a 2 + l) = ^(Ia- 3 p5 + aV(0)) (40) 

which develops a singularity corresponding to the part of the final black hole singularity 
covered by matter. 

The corresponding effective classical equations are modified by replacing the classically 
diverging a~ 3 in the matter Hamiltonian with a regular function d(a) derived from finite 
inverse scale factor operators such as (J53j) jZHllZSl- Including two ambiguity parameters j 
(a half integer) and < / < 1, this can be parameterized as 

rf(a),,:=a- 3 ^(3a 2 /7j4) 3/(2 - 20 

with 

- j^q ((? + - sgn(g - l)\q - 

The essential property of d(a)jj is that it is increasing from zero for a 2 < ^jjip which 
through the Friedmann equation implies a different dynamical behavior at small volume. 
This model then provides an intuitive explanation for the removal of classical singularities 
even at the effective level since the equations lead to a bounce at non-zero a: due to the 
modified density the kinetic term is negligible at small a, and matter evolution equations 
from the modified matter Hamiltonian imply friction of <fi [80 . The potential term is then 
dominating and almost constant which means that the bounce is approximately of de Sitter 
form [52] . In this interpretation of collapsing matter this means that it does not collapse 
completely but rebounds after a point of minimal contraction is reached. 



(41) 
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Figure 2: A closed Friedmann-Robertson- Walker model and a generalized Vaidya metric 
are matched to form a collapse model. Singularities are indicated by dashed lines and the 
matching surface by dotted lines. The bottom parts of the diagrams are cut off since they 
depend on details of the solutions. The left hand side shows the classical case with a future 
singularity, while the right hand side shows the singularity-free effective case. 

So far, we had only access to the inside of the matter contribution which we assumed to 
be isotropic. The solution can now be matched to a suitable solution describing the outside, 
which would be able to tell us, for instance, whether horizons form. For pressure-less matter 
one can match to the static Schwarzschild solution as in the Oppenheimer-Snyder model 
[HTj . This is the case classically only for dust, which however can develop pressure if 
quantum modifications come into play. (This is in agreement with our earlier observations 
that quantum effects will not allow the presence of a static vacuum solution.) We have 
thus chosen the more general scalar matter which has pressure even classically. Physically, 
pressure leads to shock waves at the outer boundary giving rise to a non-static exterior. 
This can be described by a generalized Vaidya metric 

ds 2 = -(1 - 2M(x, v)/ X )dv 2 + 2dvd X + x 2 dfi 2 (43) 

which we can match to the interior (Fig. |2J) by requiring equal induced metric and extrinsic 
curvature at the time-like matching surface S defined by r = R inside and x — x( v ) outside. 

In this way we can see what the collapsing matter implies for the outside space-time 
at least in a neighborhood 66 . To have access to the full outside all the way up to an 
asymptotic observer we would need to specify the matter content outside. While this 
is possible e.g. as the same scalar matter as inside, only inhomogeneous, it would not 
necessarily be correct physically. In fact, we have modified only the classical equations 
describing the interior, while we did not use effective equations outside. When the matter 
distribution extends over a large region in the early stages, one does not expect strong 
modifications, but this is not clear close to the bounce. In this region, the interior equations 
are strongly modified, and this is transferred to the outside via the matching conditions in a 
rather indirect way: We use the classical generalized Vaidya metric, but did not specify the 
matter content. It is effectively the energy momentum tensor which carries quantum effects 
from the interior to the outside via the matching. Prescribing the outside matter content 
would remove this transfer and stop us from seeing possible quantum effects outside. 
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We write the interior metric as 



ds 2 = -dt 2 + X(r, t) 2 dr 2 + Y(r, t) W 

with X(r,t) = a(t)/(l + r 2 /4) and Y(r, t) = rX(r,t). On the matching surface r = R 
of the interior and x = xi v ) °f the generalized Vaidya exterior the metric and extrinsic 
curvature have to agree. From the metrics we obtain 



and 



dv 
dt 



[l-2M/ X -2d X /dv)- 1 / 2 \ 



(44) 
(45) 



while the extrinsic curvature, computed again from ([21)1. gives us 



YY' 



X 



x- 



1 - 2M/ X - d X /dv 



and 



d 2 y / 2M dy 
= dv M+^4+ 1 3 A 



dv 2 



1 - 2M/ X - 2d X /dv 
M 



X 



dv 



X 



-d x M 



(46) 



(47) 



which yields a condition for dM/dx at constant v. 

With dx/df|s = xIs/^Ie and (jB]) we use (j4*o^) to write the square root in (J43|) in terms 
of V and Y which leads to 



dv 
dt 



_ (Y'/X + Y) 

s i - 2M/r 



(4E 



Using (J44j) and defining c := Y'/X, (J46|) becomes 

c 2 (l - 2M/ X + 2d X /dv) = (1 - 2M/ X - d X /di;) 2 

which with 

(dx/dv) 2 = Y 2 (l - 2M/ X - 2d X /d^) 
(following from dx/dv = x/ v an d (03) gives c 2 = 1 — 2M/x + Thus, 

2M| s = (rr 2 + r(i-c 2 ))| s . 



(49) 



A trapped surface forms in a generalized Vaidya metric when 2M = x, which lies on 
the matching surface if 2M = Y. From (J49|) this yields the simple condition 



which for FRW reduces to 



\Y\ = c = Y'/X 
= (l-R 2 /4)/R. 



(50) 
(51) 
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Figure 3: Sketch of the bouncing 
effective interior, covered by evap- 
orating horizons (dashed). Con- 
stant v (outside) and t (inside) 
slices are dotted. The matching 
described in the text only refers to 
the part before the inner horizon 
collides with the expanding mat- 
ter after the bounce. 



Assuming, for now, that this condition will be satisfied at a time t(R) during collapse, we 
obtain a horizon covering the bounce (Fig. EJ). The squared norm of its normal is given by 
d v M(l — 2d x M) which can be computed from (|4T)) using dM/dv = d v M + d x Mdx/dv and 
turns out to be zero if the horizon condition ()51|) is satisfied. The horizon is thus always 
null when it first intersects the matching surface. 

After the first trapped surface forms on the matching surface, |d| continues to increase 
before it turns around when the peak in d(a)j t i is reached. From then on, |d| decreases 
and reaches a = at the bounce. In between, the trapped surface condition (J5TJ) will be 
satisfied a second time at an inner horizon. Unlike the outer one, it lies in the modified 
regime where energy conditions are effectively violated and a > [44 . It is also null at the 
matching surface but can become time-like soon and evaporate later. Similarly, the outer 
horizon can become time-like when matter having experienced the quantum modifications 
starts to propagate through it. Thus, also the outer horizon can evaporate and shrink 
toward the matching surface at later times, when the inner matter is already expanding. 

The horizon thus evaporates and the bouncing matter has a chance to reappear. Indeed, 
the condition (J51|) will be satisfied also at a time after the bounce where a is now positive. 
Thus, the horizon will intersect again with the matter shells and disappear. At such a 
point, however, the matching breaks down since dv/dt diverges when 2M = Y and Y > 0. 
From a single matching of the interior we obtain only a part of the collapse before a horizon 
disappears. At the endpoint of the horizon the interior coordinate t ceases to be good, and 
we have to match to another patch (Fig. |HJ) • 

The precise position of the horizon only follows when we specify the matter content 
and field equations (i.e. Einstein's equations or modified ones) and integrate with M and 
dM/dx as boundary conditions at the matching surface. Since we leave this open, we 
do not get precise information on the horizon but only qualitative properties. After some 
time into the collapse, the horizon is expected to shrink since modifications in the interior 
imply small violations of energy conditions (which also allow the bounce to take place). 
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Radiation of negative energy implies, analogously to Hawking radiation, that the horizon 
becomes time-like and shrinks. Later, it can meet the matching surface again at which 
point the matter becomes visible from behind the horizon. If the initial mass was large, it 
takes a long time for the bounce to occur and the matter to reemerge such that for most 
of the time the system looks like a classical black hole to an outside observer. 

It is not guaranteed that a horizon forms at all since fulfillment of the horizon condition 
depends on initial values. In particular, once we choose R to specify the matching surface in 
the interior, Eq. (|51|) fixes the value for a which needs to be reached for a horizon to occur. 
Classically, a is unbounded as we approach the singularity such that the condition will 
always be true at one point and there is always a horizon covering the classical singularity in 
this model. With the effective equations, however, a is bounded for given initial conditions, 
and depending on the value of R it can happen that the horizon condition is never fulfilled. 
In this case, there would not be a black hole but only a matter distribution collapsing to 
high densities and rebounding. This rules out black holes of a certain type, in particular 
those of small mass: Starting with a configuration such that a horizon forms, we can 
decrease R toward zero without changing a(t). The right hand side of (|51|) then increases 
and at one point the condition can no longer be fulfilled. Since with decreasing R we carve 
out a smaller piece of the homogeneous interior, the total initial mass is smaller, giving us 
a lower bound for the mass of black holes in this model. Precise values have to be derived 
from more detailed models, but this argument shows that large, astrophysical black holes 
will be unaffected while primordial ones of small masses do not form 66J. 

5 Extrapolation 

We have seen two results from homogeneous techniques employed in the preceding section: 

• At the fully quantum level of the Kantowski-Sachs model describing the Schwarzschild 
black hole interior the singularity is absent (Fig. [TJ. 

• Matter systems allow effective classical equations for their collapse such that the 
classical singularity is replaced by a bounce sometimes shrouded by a horizon (Fig.^J. 

Both results have been arrived at with very different techniques, and have different physical 
meaning. The first one only applies to the vacuum case but provides us with a strict result 
as to how the classical singularity is replaced in quantum gravity. It directly shows that 
general relativity is singular because it relies on the smooth classical space-time picture. 
This picture breaks down at high curvature and has to be replaced by discrete quantum 
geometry, providing a non-singular evolution. 

The second result works with matter but is more intuitive, only giving a picture from 
effective classical equations. It provides a physical, rather than geometrical explanation 
for the failure of general relativity in strong curvature regimes. Singularities in general 
relativity can be understood as a consequence of the always attractive nature of classical 
gravity: Once matter collapses to a sufficiently high density, be it an isolated part or the 
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5 EXTRAPOLATION 




Figure 4: Combination of Figs. Q and 
01 where the quantum region around 
the classical singularity is shrouded by 
an evaporating horizon (dashed). The 
dotted line marks the boundary of the 
part of space-time unaffected by the 
strong curvature inside. 



whole universe, there is nothing to prevent total collapse into a singularity. Viewing the 
Friedmann equation, e.g. with scalar matter, as describing a mechanics system with the 
matter energy density serving as potential shows this by the fact that the energy density 
decreases as a function of a at fixed <fi and p^, in particular the kinetic term ^a~ 3 p^. Thus, 
there is an attractive force driving the system toward a = (or, as usually expressed in 
cosmology, positive pressure which thermodynamically is defined as the negative change 
of energy with volume). The modification of a~ 3 by the regular function d(a) in (|41j). 
which turns around at a peak value and then approaches zero rather than infinity at a = 0, 
implies that now the energy density increases as a function of volume at small scales. This 
can be interpreted as quantum gravity becoming repulsive at small scales, which can then 
easily prevent total collapse into a singularity. Moreover, at non-zero but small scales this 
repulsive component is still active and leads to modified behavior. For instance, in an 
expanding universe it implies that the expansion is accelerated leading directly to inflation 
|44j . In the interpretation of collapsing matter, the same effect makes the horizon shrink 
after the bounce such that only the strong quantum region is covered. 

Since we have used approximations, the question arises how these partial results can fit 
into a full picture of quantum black holes. The first result indicates that space-time can 
be extended through classical singularities, but since it gives us access only to the interior, 
it is not clear if the new region we reach can also be accessed from an outside observer. (If 
not, the black hole would appear as a wormhole through which one can travel into a new 
region of the universe.) The second result now indicates that we can in fact access the new 
region since there is only one matching region outside the collapsing matter, suggesting a 
picture as in Fig. 0] Quite similar ideas have been put forward, by different motivations, 

in |H21E31E3|- 

However, here it is important to bear in mind that we only effectively described matter 
outside falling into the collapsing shells. In a more realistic model, there would be such 
inhomogeneous matter colliding with the homogeneous core and making it more heavy. 
This can then lead to singularities forming in the outside region. For their resolution we 



25 



would again have to use quantum geometry and face the same problem as to whether or 
not this will lead to a new region accessible from the outside. 

It is clear that a decisive answer can only be obtained with inhomogeneous techniques, 
which we are going to describe in the next section. Still, even at this level one can see that 
there are only a few possible scenarios which can be distinguished by using inhomogeneous 
properties of quantum geometry. Irrespective of which outcome inhomogeneous models 
will show us, one can already see special features of quantum geometry leading to a new 
paradigm about black hole evaporation. For the first time, this takes into account a 
resolution of the classical singularity with implications for apparent loss of information 
In fact, while Hawking radiation still emerges in a neighborhood of the dynamical 
horizon and is still approximately thermal, this is by no means everything coming out of 
the black hole at late times. Infalling matter now evolves through the quantum region 
of high curvature and reappears later, restoring correlations which are not recovered by 
Hawking radiation alone. In particular, there is no reason for the final state measured on 
all of future null infinity to be mixed if we started with a pure initial state. In the usual 
picture one would cut out the quantum region (or the place of the classical singularity) and 
consider the future space-time without allowing penetration through that region. Future 
null infinity then stops at the intersection with the dotted line in Fig. 01 and a state 
retrieved at this part of null infinity is indeed mixed since it is obtained by averaging over 
the rest to the future. In this way, both the singularity problem and the information loss 
paradox are resolved by loop quantum gravity. 

6 Inhomogeneous techniques 

For the spherically symmetric model we need to perform the loop quantization for inhomo- 
geneous configurations (|17|) and (jllj) such that the basic fields now depend on the radial 
coordinate x. Instead of using states such as ()27|) with a finite number of labels, we now 
have a field theory with infinitely many kinematical degrees of freedom. An orthonormal 
basis of states is given by [T3] 

(A X ,K V \ ...,k n ,n n , k n+1 , fj, n+1 , . . .) = J|exp(|iA; n / A x dx) ex^{-ip, n ^K v (v n )) (52) 

n 

with countably many labels k n G Z and < fi n G K labeling edges e n and vertices v n , 
respectively, of a 1-dimensional graph in the radial line. Note that, as already indicated 
before, we are using exponentials of the extrinsic curvature component along homoge- 
neous directions but holonomies of the connection component A x along the inhomogeneous 
direction. Both exponentials are represented as multiplication operators. 
Spatial geometry is encoded in densitized triad operators acting by 



E*(x)\ 
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where n (x) is the edge label to the right (left) of x, and X is an interval on the radial line 
(over which we need to integrate E v since it is a density). As before in the homogeneous 
case we also obtain densely defined operators for inverse powers of the triad components, 
which can in particular be done for the inverse of E v in the spin connection component 

6.1 Hamiltonian constraint 

For the Hamiltonian constraint (|16|) we again have to obtain curvature components from 
holonomies, now using holonomies of A x for the inhomogeneous radial direction and ex- 
ponentials of K v for homogeneous directions along symmetry orbits. Terms containing 
the spin connection component T v belonging to homogeneous directions will be quantized 
separately. One may wonder if this procedure will easily give the right components in 
the Hamiltonian constraint, given that it has the rather simple expression fll6j) in terms 
of extrinsic curvature components while we are using the Ashtekar connection component 
A x . As we will see, the general scheme will yield automatically the right combination of 
components by a straightforward construction of loops to be used in holonomies. 

To see this in detail, we first note the difference between A x and K x , which is given by 
the x-component of the spin connection. For a general spherically symmetric triad it takes 
the form 

E x/ E xl 

T = — r?V 3 dx H Ad$ A sin -d&ip + cos ^do? (55) 

2E* 2E Lp y ' 

as in Q with (fTU|) . Here, we recognize (|T^|) as used before as the component along homo- 
geneous directions. This component is a scalar (noting that both E x ' and E v are densities 
of weight one), while the x-component T x = —rf does not have covariant meaning and 
indeed can be made arbitrarily small locally by a suitable gauge transformation. This 
is analogous to the situation in the full theory, while the gauge invariant meaning of 
mimics homogeneous models. 

Even though r x can be made arbitrarily small locally by a gauge transformation, we 
cannot assume this when constructing a suitable Hamiltonian constraint operator. Thus, 
it must be built into the construction so as to combine with A x from radial holonomies 
to give K x as in the expression for the constraint. This K V K ^-term in the constraint can, 
according to the general construction where connection or extrinsic curvature components 
derive from closed holonomies, only come from a loop which has one edge along a symmetry 
orbit and one in the radial direction. Starting in a point v, such a holonomy is of the form 
h x ^h^\v + )(h x S ^)~ l (hy\v))~ l with a new vertex v + displaced from v by a coordinate 
distance S x of the radial edge. (We distinguish between S x for the radial direction and 5 
for the angular directions since the continuum limit is technically different in both cases.) 
This term appears together with A(v) (coming from quantizing the triad components) in 
a trace whose expansion in 5 

-2ti(h x h^v + )h- l h^v)- l K(v)) 
~ ^(K^tTiAiv^Aiv)) + 2K v (v+) J A x dxtT{r 3 A(v+)A(v))) 
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= j8(K ip (v+) sin(?7 + -T]) + K ip (v+) J A x dx cos(?7 + - n)) 

= 88 xl K v {v){A x {v) + rj'(v)) + 0(5 2 ) (56) 

has all the right terms, with A x coming directly from the radial holonomy and i](v + ) in 
5i] f ~ i]{ v +) ~ v( v ) from the internal direction A(v + ) at the new vertex. The other term of 
the form K% is obtained from angular holonomies only, as in the homogeneous case, 

- 2 tr(^^^ 1 / i - 1 r 3 ) ~ 5 V^J • (57) 

The matrices T3 and A(v) in the traces come from Poisson brackets expressing triad com- 
ponents as in the homogeneous constraint. Moreover, the spin connection components in 
(fT6*|) are again expressed through the curvature 

F(T) = -T'^Adx Adi? + r^Asinte A dip + (rj - l)r 3 sin« A dip 

of the spin connection. Through 

t IJK tr^jF^jjhKih-K 1 , V}) ex (Tj - 1){A X , V} - 2r;{- 7 ^, V}) 

we obtain the additional terms of the constraint, where we included the length parameters 
5 1 (i.e. 5 X for I = x and 5 for $ or <p). 

This demonstrates how the general procedure works without additional input: We use 
exponentials of extrinsic curvature components for homogeneous directions and holonomies 
of Ashtekar connection components for inhomogeneous ones as dictated by the background 
independent representation. Spin connection components for inhomogeneous directions 
then come in the right form to combine with extrinsic curvature components, while those 
in homogeneous directions are split off and quantized separately. This is possible because 
those components, in contrast to the inhomogeneous ones, do have covariant meaning. This 
ties together the constructions in homogeneous models and the full theory, and at the same 
time opens a direct route to effective classical equations: Homogeneous spin connection 
components usually contain inverse powers of the densitized triad, such as (|13p. When 
they are quantized, the classical divergence will be removed implying modifications at 
small scales, in homogeneous models this has been used, e.g., in the Bianchi IX case 
where it has been shown to remove the classical chaos [H3 [HH] - Similarly, one can use this 
mechanism to derive effective classical equations for the spherically symmetric model and 
find possible consequences. 

Before trusting those effective equations in the inhomogeneous case one needs to make 
sure that there is a well-defined Hamiltonian constraint operator emerging from the proce- 
dure described here. So far, we have only discussed those holonomies and spin connection 
components which give us the contributions to the constraint, but they must now be stuck 
together with quantizations of triad components so as to build a well-defined operator for 
the whole expression. Moreover, since the expression is an integrated density, one has 
to discretize the integration first and then, after quantizing the individual terms, perform 
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the continuum limit removing the regulator. The discretization had already been under- 
stood above, with v and v + being the endpoints of a discrete interval of size 8 X , and we 
convinced ourselves that the continuum limit of the discretization will yield the correct 
result. In more detail, one writes 

H[N] = j dxN(x)H(x) ~ ^28V>N(v n )H(v n ) 

where we discretized the radial line into intervals of coordinate length 8 X , each one con- 
taining the point v n . Classically, both expressions can be made to agree for any subdivision 
by choosing points v n in the intervals according to the mid point theorem. Alternatively, 
if one wants to fix the v n to be endpoints of the intervals, the discretization agrees with 
the classical constraint in the continuum limit in which n — > oo and 8^ — > for all n. 

For the Hamiltonian constraint in general each term in the sum then has contributions 
of the form 

8 x H(v) oc e IJK ti(huh K [h K \V}) (58) 

I,J,K 

where we sum over triples (J, J, K) of independent directions which in symmetric models 
are given by generators of the symmetry transformations (($, if) in the spherically symmet- 
ric case), and in the full theory or inhomogeneous directions of a symmetric model by edges 
of a graph (x in spherical symmetry). The holonomies hu are formed according to the 
symmetry: if both directions / and J are inhomogeneous, hu is a holonomy along a closed 
loop ajj constructed from edges in the / J-plane of a graph to act on; if at least one of the 
two directions is homogeneous, hjj = hjhjhj 1 ^ 1 — j 2 8j8jFu with hj being exponentials 
of the (su(2)-valued) extrinsic curvature components belonging to the /-direction for a ho- 
mogeneous direction I or holonomies along an inhomogeneous direction. (The appearance 
of F(T) can be understood as a correction term since loops made from holonomies along 
vector fields generating symmetries do not close if orbits have non-zero curvature [8*7].) 
The size of loops ajj or single holonomies is determined by the size 8 X of the discretiza- 
tion. The final holonomy hx either belongs to an edge transversal to both directions / and 
J, or again to exponentiated extrinsic curvature components if K is homogeneous. These 
combinations are chosen in such a way that hu yields the correct curvature components, 
the commutator gives the necessary triad components, and both terms together provide 
just the right product of lattice sizes such as 8 X in order for the sum to take the form 
of a Riemann summation of the original integral. In the full theory, this procedure only 
results in the so-called Euclidean part of the constraint which can be used to construct 
the Lorentzian constraint [71]. In the models used here, however, the prescription (J58j) is 
sufficient even for Lorentzian signature. 

This can now be illustrated and applied in the spherically symmetric model where the 
first case above cannot appear since there is only one inhomogeneous direction. We have 
thus two cases, one in which direction / or J is radial, resulting in the first product of 
holonomies discussed above, combined with a commutator 

h^ [h^ 1 , V] = V — cos IjKp V cos ^jK^ - sin ^jK^V sin ^jK v 
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— 2A(cos rfyKtpV sin §7!^ — sin ~jK v V cos §7!^) • 

In the second case we have the other product of holonomies and a commutator 

h x [h-\V] = V-cos±fA x VcoslfA x -smifA x VsmiJA x 

+2r 3 (cos \ / sin \j A x - sin \ J A X V cos ± / A x ) . 

The integration here is over an interval of size 5 X such that in the limit of a fine discretization 
this term is of order S x as needed for the Riemann sum. 

What we did not specify yet is how the discretization is adapted to a graph the con- 
structed operator is supposed to act on, i.e. whether v and maybe v± are already vertices 
of the graph or arbitrary points. At this point, choices need to be made which lead to 
different versions of the constraint. The same choices arise in the full theory |HH1 IH9*1 19*U]. 
but they can be studied much more easily in the spherically symmetric model such that it 
may be possible to rule out some versions. 

It is already non-trivial to check that different versions lead to well-defined operators 
at all. For this, the action after performing the continuum limit, in which the number 
of discretization points becomes infinite, must be finite. If the action were non-zero at 
each discretization point, there would not be a well-defined operator in the limit and the 
regulator could not be removed. One would then only deal with a lattice regulated theory 
rather than a quantization of the continuum theory. In the full theory, there is a well- 
defined operator because the action of the constraint is zero unless a discretization point 
is already a vertex. Starting with states with finitely many vertices then leads to a densely 
defined operator. This comes about in the full theory because the constraint contains the 
volume operator in such a way that it acts only on planar vertices if there is no vertex 
already present in the graph. Since the full volume operator annihilates all planar vertices, 
there are only finitely many contributions from the vertices already present. 

In spherical symmetry, however, all vertices are planar since graphs are just 1-dimensional. 
This simple general argument is thus not available and it is not obvious that the same con- 
struction scheme will result in a well-defined operator. It turns out, however, that this 
is the case as a consequence of how triad components in the constraint are quantized: a 
discretization point which is not already a vertex of the graph to act on will be annihilated 
such that only finitely many contributions from the vertices remain. One can thus use the 
same type of operator, just with adaptations to the symmetric situation. 

Nevertheless, one can also choose different constructions where the discretization is 
given directly by the graph, i.e. discretization intervals would be complete edges of the 
graph. Since endpoints of discretization intervals are always vertices of a state, the con- 
tinuum limit would then require also states to change and become finer and finer. In this 
picture, the continuum limit of the constraint operator can only be tested on states which 
are suited to the continuum behavior, while there are also other states where discreteness 
is essential and where the classical constraint would be corrected from quantum effects. 
Moreover, in the continuum limit the number of vertices diverges and the constraint opera- 
tor becomes ill-defined just as the usual Wheeler-DeWitt operator is. Both schemes result 
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in well-defined operators, but they lead to quite different equations of motion and require 
different conceptual viewpoints about the continuum limit. When the continuum limit is 
to be ensured for each state, one requires in a sense that the classical equations are sensible 
at arbitrarily small scales, and corrections could only come from quantum uncertainties. 
In the second picture, on the other hand, the classical continuum picture arises only after 
a certain coarse graining, or by working only with states which are not sensitive to the 
microscopic structure. If one chooses a state which is sensitive to small scales, or looks 
very closely at small scales of even a semiclassical state, then corrections to the classical 
expressions arise, for instance as a consequence of the underlying discreteness. This second 
viewpoint has been taken successfully in cosmological models, and is, as we will see, also 
fruitful in black hole models. 

6.2 Dynamics 

Since the constraint operator is again constructed from holonomies which act by shifting 
the labels, it implies difference equations for states in a triad representation. (Note that 
also in the spherically symmetric model the triad representation exists, unlike in the full 
theory where flux operators do not commute These equations are now not only 

partial difference equations but also have many independent variables. Interestingly, the 
type of difference equation is very different for the two versions of the constraint operator: 
in the second case the number of edges and vertices of the original state is unchanged and 
the operator only acts on the labels. This results in difference equations with independent 
variables k e for each edge and /i„ for each vertex. Since the operator does not change the 
number of edges and vertices, one obtains coupled difference equations in a fixed number 
of variables for each sector given by the number of vertices. 

In the first version of the constraint, however, the situation is very different. Now, 
new vertices are created and edges split in each action of the operator. Thus, for a triad 
representation it is not enough to work with a fixed number of vertices. Rather, all graphs 
have to be taken into account for the equation, which implies that one has to deal with 
infinitely many independent variables and thus functional difference equations. 

We thus return to the simpler type of difference equation implied by the other version 
of the constraint and discuss what one can already say about the singularity issue. First of 
all, one will have to identify classical singularities on minisuperspace in order to study the 
constraint equation in a neighborhood. In an isotropic model this is simple since the only 
way is for the volume to go to zero ^Uj- Similarly, one can identify the classical singularity 
on an anisotropic minisuperspace where all densitized triad components would go to zero. 
The situation is not so clear in midisuperspaces such as spherical symmetry since there are 
more possibilities for a singularity to develop. Even though this has not been settled in 
general, there are many cases where a singularity is characterized by E x approaching zero 
(which on classical solutions such as Schwarzschild implies that becomes zero, too). 
This is also in agreement with the general mechanism removing singularities seen so far in 
homogeneous models: there it is the sign coming from orientation which leads to different 
regions of minisuperspace separated by the classical singularity. The quantum evolution, 
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as we have seen in the Kantowski-Sachs model, can then allow us to evolve between the 
two regions, thus removing the classical singularity as a boundary. 

The role of orientation is now played by sgnE x since det q = E x {E ip ) 2 . Since E x 
depends on the radial position, or the edge after quantization, the boundary of our midis- 
uperspace has many components which we can identify with an inhomogeneous classical 
singularity, corresponding to the fact that inhomogeneous singularities behave differently 
in their different points. For the states this means that we encounter a section of a classical 
singularity each time an edge label k e becomes zero. As in homogeneous models, we can 
then use the evolution equation in the triad representation in order to see if an evolution 
through this part of the boundary is possible. Since the structure of the difference equa- 
tion for a given edge label is very similar to the homogeneous equation, one can expect 
that the boundary indeed disappears and that the quantum evolution connects regions 
of midisuperspace corresponding to different local orientations. There would thus be no 
singular boundary, and the same mechanism as in homogeneous models could also remove 
spherically symmetric singularities. 

This scenario has been verified in noting a crucial difference to homogeneous mod- 
els wich require a symmetric ordering of the constraint. Thus, quantization choices are 
reduced by looking at less symmetric models, so far in such a way which maintains the 
validity of the general picture. Yet, there are also open issues left for a general under- 
standing. For instance, while the results are independent of the matter Hamiltonian and 
can be extended to cylindrical gravitational wave models, thus also allowing local degrees 
of freedom, they are so far based only one type of the constraint which leads to difference 
equations easier to deal with. The behavior with the other version is not easy to see, but 
if it is singularity free, too, the mechanism is likely to be different. Most importantly, the 
kind of initial/boundary value problem suitable for the constraint equations needs to be 
analyzed in more detail to guarantee that there are suitable and sufficiently many solutions 
with the correct classical limit. At this point the anomaly issue, i.e. whether two constraint 
operators with different lapse functions have the correct commutator, becomes important. 
These considerations thus provide a promising and treatable way to distinguish different 
versions of the quantization by their physical implications, which can then be extrapolated 
to the full theory. 

6.3 Horizons 

A feature of black holes which is new compared to cosmological models, and which requires 
inhomogeneous situations, is given by the presence of horizons. Global concepts such as 
the event horizon are, of course, not helpful in our case since we would need to solve the 
Hamiltonian constraint completely before being able to discuss this issue. There are more 
practical definitions such as apparent horizons which, however, are much more general 
and do not distinguish between fully dynamical situations and almost static systems. The 
quantum behavior would be most easy to analyze if we can define horizons locally and 
in a controlled manner which does not require the full dynamics at once. Such a concept 
is presented by isolated |2j or slowly evolving dynamical horizons which even quite 
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unexpectedly simplify the spherically symmetric Hamiltonian constraint in their neighbor- 
hood. 

6.3.1 Definition 

There are three main parts to the definition of an isolated horizon A with spatial sections 
S = S 2 of given area do, embedded into the space manifold E by i: S — ► £ ITTH] : 

(i) The canonical fields (A l a , E?) on the horizon are completely described by a single 
field W = \i* A % ti on S which is a U(l)-connection obtained from the pull-back of 
the Ashtekar connection to S. Here, is an internal direction on the horizon chosen 
such that W is a connection in the spin bundle on S 2 and r l E°; = y/det q r a on the 
horizon with the internal metric q on S and the outward normal r a to S in E. 

(ii) The intrinsic horizon geometry, given by the pull-back of the 2-form := e a b c E^ to 
S, is determined by the curvature F = dW of W by 



(iii) The constraints hold on S. 

A further consequence of the isolated horizon conditions jHS] is that the curvature T of 
the pull-back of A\ to S has to equal the curvature of W: rjjF(i* A 1 ) = 2dW. This can be 
seen as one of the distinguishing features of an isolated horizon since even slowly evolving 
horizons at rate e (related to the expansion of horizon cross sections |92j ) will break it, 
though just by an amount of the order e 2 . 

When the horizon is introduced as a boundary, condition (i) is used to identify the 
horizon degrees of freedom represented by the field W. Condition (ii) then shows that 
these degrees of freedom are fields of a Chern-Simons theory on the horizon. It is the 
main condition since it relates the horizon degrees of freedom to the bulk geometry, which 
after quantization selects the relevant quantum states to be counted. Condition (iii), on 
the other hand, does not play a big role since an isolated horizon as boundary implies a 
vanishing lapse function on S for the Hamiltonian constraint which then is to be imposed 
only in the bulk. 

Thus, when computing black hole entropy in this way, as we will describe later, the 
Hamiltonian constraint does not play any role since it does not act at the boundary, and 
the Hamiltonian generating evolution along the horizon need not be considered. In the 
spherically symmetric model one can hope that the constraint is simple enough for an 
application in this case, either to generate evolution or to impose the horizon not as a 
boundary but inside space such that the constraint would have to be imposed. In the 
latter case, moreover, we will not be able to have an independent boundary theory which 
is then matched to the bulk, but would have to find the relevant degrees of freedom within 
the original quantum theory. 



F 




(59) 
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6.3.2 Spherical symmetry 

We can now evaluate the conditions for spherically symmetric connections of the form 

A = A x (x)r 3 dx + A (p (x)A A (x)d'd + A (p (x)A A (x) shindy? + r 3 cos$oV (60) 

and densitized triads (jllj) . where in general the internal directions A A and A are different. 
The connection component A x has been discussed before, while the relation A^A A = 
T^A — '-fK^A, following from the definition of the Ashtekar connection together with Q, 
(JTOJ) and flni), implies A 2 ip = Y 2 + ^K*. 

We choose r« := sgn(E x )Si^ such that in fact r % E^ = \E X \ siwddx with the intrinsic hori- 
zon area element l-E^I sin^ of a metric jE^ldfi 2 . Thus, W = \rii*A % = | sgn(E x ) cosfMy? 
whose integrated curvature given by <f s dW = —2n sgn(E x (x)) agrees with the Chern num- 
ber of the spin bundle, depending on the orientation given by sgn(E x ). 

Evaluating (J59|) first shows that in the spherically symmetric context it is not restrictive 
since we have a = 4ir\E x (S)\ and the right hand side given by — ~ sgn(E x (S)) equals F 
for all E. This is not surprising since the spherically symmetric intrinsic geometry of S 
is already given by the total area which is fixed from the outset. (What is free is the 
sign of E X (S), or orientation, which confirms ideas of [94J.) Now the first condition plays 
a major role, which we evaluate in the form riF^CA 1 ) = 2dW jHEj- Since ^{CA) = 
(A 2 — 1)t3 sin^d'i? A dip, the condition requires A v = which will be the main restriction 
we have to impose on quantum states in addition to the constraints. This condition A^ = 
selects 2-spheres in a spherically symmetric space-time corresponding to cross-sections of 
a horizon. Indeed, for the Schwarzschild solution we have A^ = since the extrinsic 
curvature vanishes. With ()13|) and the Schwarzschild triad we obtain the correct condition 
x = 2M for the horizon. In general, A 2 Lp = Y 2 ip + ^y 2 K 2 = implies = and K v = 0. 

A slowly evolving horizon at rate e satisfies the condition rjjF(t*y4*) = 2dW only up to 
terms of the order e 2 . Thus, A v is not exactly zero but must be small of order e, which 
then is true also for and K^. 

6.3.3 Dynamics 

In spherical symmetry we can locate a horizon on a state [69J, which must be at a vertex 
in order for = to be sharp enough. The condition that K v be zero for an isolated 
horizon or small for a slowly evolving horizon then leads to important simplifications which 
allow a perturbative treatment of the dynamics around the horizon. Indeed, when acting 
with the constraint at the horizon both terms made from holonomies hjj contain factors 
of sin at the horizon vertex or a neighboring one which must then be small. Ignoring 

those terms in an approximation leads to an operator which is diagonal on the spin network 
states ()52|) and thus easy to solve as a constraint or to use for generating time evolution 
at a boundary. The additional terms ignored in this approximation can then be included 
in a perturbative treatment of the near horizon dynamics. 

Without many calculations this already shows how the horizon fluctuates dynamically. 
Classically an isolated horizon has constant area which thus commutes with the Hamil- 
tonian constraint. This is also true at the quantum level to leading order of the above 
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approximation since the area operator A(S) = Att\E x (S)\ has the same eigenstates (|52|) 
as the leading order constraint. Thus, at this level the horizon area is an observable not 
just when the horizon is treated as a boundary, but also if its full neighborhood is quan- 
tized. However, there are additional terms which arise in higher orders of the perturbation 
scheme. There are two reasons for horizon area fluctuations even in the isolated case: While 
classically K v = exactly at the horizon and only this value is important, the quantization 
does not allow this to hold arbitrarily sharply. Otherwise, the volume of a shell around 
the horizon, which depends on the conjugate momentum E v of K v could not be sharp 
independently of the mass which would contradict semiclassical properties to hold true at 
least for massive black holes. Secondly, the constraint operator acting at the horizon itself 
depends on neighboring values of K v through h v (v + ) in (fo7)|) . This would give non-zero 
contributions even if K v at the horizon would be zero exactly. 

Both terms lead to small dynamical changes in the horizon area coming from typical 
quantum gravity properties. The first reason is quantum uncertainty which does not allow 
a sharp condition = 0, and the second space-discreteness and non- locality which implies 
that not only K v at the horizon itself is relevant but also the values in neighboring vertices 
which are not necessarily zero. For large black holes, the correction terms are expected 
to be small: uncertainty will not change the horizon area much compared to its already 
large size, and in neighboring vertices of a semiclassical state K 9 will still be extremely 
small. Thus, for large black holes the horizon area is an excellent observable, while for 
microscopic black holes large fluctuations are expected which may even prevent horizons 
as they are known classically. This agrees with the picture we have obtained from effective 
equations and matching techniques before. 

7 Full theory 

The methods developed so far in symmetric models mimic those of the full theory, with 
some adaptations to preserve the symmetry. In this section, for completeness, we describe 
what this looks like in the full theory and discuss applications which work without assuming 
symmetries. 

7.1 Representation 

As discussed before, the full theory of loop quantum gravity is based on holonomies for 
arbitrary edges in space and fluxes for surfaces, forming the basic classical Poisson algebra. 
In the connection representation, states are functionals on the infinite dimensional space 
of connections jHSl ESI ED] through holonomies, and a dense subspace of the Hilbert space 
is spanned by cylindrical functions 

1>(A) = Mh eL {A),...,h en (A)) (61) 

which depend on only finitely many holonomies. Since there is now no symmetry require- 
ment, the edges can be arbitrary curves in space and form a graph 7 with vertices at their 
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intersection points. The inner product for two states associated with the same graph is 
given by 



with the Haar measure dfin on the structure group SU(2). For two functions with different 
graphs, they need to be extended to a bigger one which is always possible by cutting edges 
or inserting new ones on which the extended state depends trivially An orthonormal basis 
is given by spin network states [H3IIIE], associated with graphs labeled by irreducible SU(2) 
representations j e at edges and contraction matrices C v at vertices, of the form 



where the representation matrices pj e (h e (A)) evaluated in edge holonomies are multiplied 
together in vertices according to the symbols C v . 

Fluxes are quantized as derivative operators in the connection representation since the 
densitized triad is conjugate to the Ashtekar connection. Replacing the triad components 
in (JT5|l by functional derivatives and acting on a cylindrical function, we obtain 



which has contributions only from intersection points y of the surface S of the flux with the 
graph 7 associated with the state. Moreover, each derivative operator for an intersection 
point can be seen to be equivalent to an angular momentum operator such that its spectrum 
is discrete and equidistant. Since there is a finite sum over all such contributions, the 
spectrum of flux operators is discrete, too. Not all the angular momentum operators 
involved necessarily commute, and so triad operators do not always commute with each 
other such that a triad representation does not exist j^I] (unlike in the symmetric models 
studied before). 

The densitized triad describes spatial geometry, and spatial quantum geometry is en- 
coded in flux operators. From the basic ones one can construct geometrical operators such 
as the area jH2HHB| or volume operator [31] which also have discrete spectra. Thus, quantum 
spatial geometry is discrete in a precise way, given by the spectra of geometric operators. 
The area spectrum is known completely, but for the volume operator this is impossible to 
compute explicitly since arbitrarily large matrices would have to be diagonalized. Spatial 
geometry at the quantum level is thus rather complicated in general if explicit calculations 
need to be done. 

This translates to the Hamiltonian constraint and other operators, for which the volume 
operator plays a crucial role. Classically, the Hamiltonian constraint is given by j2Ij 




(62) 



T^,c(A) = l[C v .l[p je (h e (A)) 



(63) 




H[N] = (811G)- 1 / d 3 xA^(x)|det 



E 



2 (K b E«E b k e ljk - 2(1 + l^KlK^EZE*) (64) 
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with the curvature F l ah of the Ashtekar connection, and the extrinsic curvature K l a = 
7~ 1 (Tj l — AD a function of the basic variables through ©• Both parts of the constraint 
can be quantized using building blocks similar to (}58|) . resulting in a well-defined operator 
[7T] even when matter Hamiltonians are included jHZj- Edges for the holonomies have to 
be chosen, which can be done in a diffeomorphism invariant manner and even in such a 
way that the quantization is anomaly free at least on states satisfying the diffeomorphism 
constraint [55] . 

Which version of the quantization is the correct one, however, is still an open issue 
since in particular the classical limit and that of perturbations on a classical background 
( "gravitons" ) are difficult to analyze. Moreover, finding and interpreting solutions in full 
generality is complicated by technical and conceptual problems. 

It is thus important to devise approximation schemes, other than symmetry reduc- 
tion as employed before, in order to shed light on physical properties of the full theory. 
One powerful possibility consists in imposing an isolated horizon as a boundary since 
boundary conditions imply that the constraint is not to be imposed there (a constraint 
has lapse function going to zero at the isolated horizon). Thus, also at the quantum level 
the constraint operator can be ignored and aspects of the basic quantum representation 
receive physical meaning. Indeed, boundary degrees of freedom are obtained from intersec- 
tions of spin network states with the horizon surface, and flux operators are important to 
select physical states corresponding to an isolated horizon. By counting those states and 
comparing with the Bekenstein-Hawking expectation the theory can be tested. 

7.2 Black hole entropy 

An isolated horizon S with prescribed area do as a boundary leads to an additional bound- 
ary term in the symplectic structure p^. 

tt = (87T7G)- 1 / d 3 xdA[ A dE? + — (167T7G)- 1 I dVi»j0A* A dA{e ab (65) 

J 27T J S 

where we denote differentials on field space by d, e ab is the anti-symmetric tensor on 
the boundary surface, and rj the internal vector as in the definition of an isolated horizon. 
The boundary term to the symplectic structure can be recognized as that of a U(l) Chern- 
Simons theory which thus describes the horizon degrees of freedom by the U(l) connection 

W a = \nAl 

We quantize the full system by using quantum geometry in the bulk and quantum 
Chern-Simons theory on the horizon. Doing this results in the curvature F = dW be- 
coming an operator with equidistant spectrum which, via (JSTJj) needs to be matched to the 
flux rjS* through the horizon. As shown before, quantum geometry indeed implies a flux 
operator with equidistant spectrum such that the matching is possible at the quantum 
level. Since also the pre-factors match, there are always solutions to the horizon condition 
which can now be counted, for a given area a , to compute the entropy as the logarithm 
of the number of states. 
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This results in an expression for entropy which is proportional to the horizon area [7J |H] , 
confirming expectations from semiclassical considerations. Intuitively, entropy counts the 
number of ways that one can construct a macroscopic horizon of area ao from elemen- 
tary discrete parts [100 (which is generalized in this picture since there is not just one 
elementary type but different ones given by the spin label of an intersection point with 
a spin network). Since the discreteness scale is set by the Barbero-Immirzi parameter 7, 
the number of possible such configurations and thus entropy must depend on 7. Indeed, 7 
appears in the constant of proportionality between entropy and area which allows us to fix 

7 by requiring the Bekenstein-Hawking law. Moreover, since there are different types of 
black holes — charged, distorted, rotating or with non-standard matter couplings — and 
the value is already fixed by the simplest case of a Schwarzschild black hole, one can test 
the theory since now entropy must result in the right way without any further parameter 
to tune. This is indeed the case |101| 1102] . providing a non-trivial test of the theory. 

The scale of discreteness is then fixed which, since it must be small enough, can already 
be confronted with observations. It turns out that 7 = 0.2735 103. 1104j is of the order 
one such that the discreteness lies around ~ \i? and is thus much too small to be 

observable directly. Indeed there have long been reasons to expect a scale of discreteness 
around the Planck length which is now confirmed by detailed calculations in loop quantum 
gravity. It is not at all obvious that this comes about since there are many non-trivial 
steps in the derivation, and mistakes in the foundations of the theory could easily lead to 
larger values which could already be in conflict with observations. 

At this point it is important to consider the physical meaning of 7. It can be seen as 
a fundamental parameter setting the scale of discreteness which is thus characteristic of 
quantum gravity. (In fact, one can express the continuum limit as a limit 7^0 [77].) 
In usual arguments, this is expected to be done by £p, which has to appear anyway just 
for dimensional reasons. However, in £p only the gravitational constant G and Planck's 
constant h enter such that the Planck length is already fixed by classical gravity and 
quantum mechanics alone. Since these theories are unrelated to full quantum gravity, 
there is no reason for 7 to equal one even though one can expect a value of the order 
one from dimensional arguments. A precise value for the scale of discreteness can only 
come from a detailed quantum theory of gravity and calculations which are sensitive to 
the underlying discrete structure, as realized by loop quantum gravity. 

8 Conclusion 

In the preceding sections we described the current status of what black holes look like from 
the viewpoint of non-perturbative, background independent quantum gravity. There are 
results obtained with different approximations to the full theory which provide a consistent 
picture of black holes without pathologies or puzzles, such as the singularity problem 
or the information loss paradox, perceived earlier from general relativity alone or from 
combinations of classical gravity and quantum field theory on a background. 

The main type of approximation used here is that of a symmetry reduction as often 
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employed in classical or quantum physics. This allows to study the background independent 
quantum dynamics and its characteristic features in different explicit ways. Compared 
to the full theory, there are several technical simplifications for instance from a volume 
operator with explicitly known spectrum. But also at a conceptual level, the interpretation 
of solutions or physical situations is simplified. 

Even though special properties of a given symmetric model, such as simplifying coor- 
dinate or field transformations, have not been made use of and essential ingredients have 
rather been modeled on the full theory, the question arises what one could do without sym- 
metry assumptions. For a fair judgment one has to bear in mind that background indepen- 
dence in quantum field theory is a new concept, which is introduced non-perturbatively. 
There are hardly any comparable results in other realistic quantum field theories, and 
quantum gravity introduces its own conceptual issues to the theory. Moreover, the fact 
that common perturbative approximation schemes are not available is a consequence of the 
property of gravity that a split into a free field theory plus perturbations is not possible. 
One thus has to deal with the fully non-linear framework which otherwise is usually avoided 
in quantum field theory. Loop quantum gravity provides a framework in which these hard 
questions, which sooner or later will have to be faced by any approach to quantum gravity, 
are being confronted directly. 

A consequence of the non-linearity is that operators, even if they can be defined in 
a well-defined manner, are by no means unique since there are often factor ordering or 
regularization choices. Loop quantum gravity, nevertheless, has succeeded in finding char- 
acteristic effects from a background independent quantization. Details, of course, depend 
on several quantization choices, but one can directly investigate the robustness of results to 
ambiguities. As described here, this allows one to solve conceptual problems in the physics 
of black holes, and also in cosmology as detailed elsewhere |1(J5| I7§1 HUt>j , while parameters 
can be fixed in detail by consistency conditions or phenomenology. 
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